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Abstract 

We calculate the running of low-energy neutrino parameters from the bottom up, 
parameterizing the unknown seesaw parameters in terms of the dominance matrix R. 
We find significant running only if the R matrix is non-trivial and the light-neutrino 
masses are moderately degenerate. If the light-neutrino masses are very hierarchical, 
the quark-lepton complementarity relation 8 C + #12 = vr/4 is quite stable, but $13,23 
may run beyond their likely future experimental errors. The running of the oscillation 
phase 5 is enhanced by the smallness of #13, and jumps in the mixing angles occur in 
cases where the light-neutrino mass eigenstates cross. 
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1 Introduction 



Low-energy neutrino experiments |]Q are providing crucial insight into lepton masses and 
mixing, though this is still limited in its scope. The most economical model for light 
neutrino masses is the seesaw model [2., but even the minimal model with three heavy 
singlet neutrinos contains a total of 18 parameters in the neutrino sector [HI E|. So far, 
neutrino experiments provide us with measurements of only four of these 0: two squared- 
mass differences Amf 2 , Am^, and two neutrino mixing angles #23,12- There are prospects 
for measuring one more mixing angle, #13 and the CP-violating Maki-Nakagawa-Sakata 
phase P 5, as well, perhaps, as the overall neutrino mass scale in cosmological data [JJ and 
one combination of Majorana mass parameters in neutrinoless double-/? decay [Bj. However, 
even these measurements would fall short of providing complete information on the full set 
of nine parameters that are in principle observable in low-energy neutrino experiments |HJ, 
out of the full total of 18. 

Nevertheless, the information available from low-energy neutrino experiments is already 
striking The atmospheric mixing angle #23 is close to maximal: sin 2 2#23 = 1.02 ± 
0.04, and the solar mixing angle 9u is quite large: tan 2 6*12 = 0.45 ± 0.05. It therefore 
seems that neutrino mixing must be qualitatively different from the smaller mixing visible 
between the left-handed quarks, where the largest mixing angle is the original Cabibbo 
angle: sin8c = 0.22. Such a difference in the quark and neutrino mixing patterns was not 
widely expected before the experiments, and has given rise to much theoretical discussion 
and speculation [TQJ dU H2] • 

One of the problematic issues in the interpretation of the low-energy neutrino data is the 
running of neutrino masses and mixing parameters below and above the seesaw mass scales. 
This renormalization inevitably introduces some 'fuzziness' in the comparison between low- 
energy measurements and any specific Ansatz for the mass matrix at the seesaw scale, since 
the renormalization depends on many of the unknown parameters in the seesaw model. The 
renormalization group equations (RGEs) have been used to study this running extensively, 
both numerically and analytically 1 . As a result, the observed low-energy neutrino mixing 
can be obtained starting from either a bimaximal Jl] or from an almost diagonal ^3] 
neutrino mass matrix at the Grand Unification (GUT) scale Mqut- I n such a situation, 
understanding the systematical features of the running of neutrino parameters becomes 
crucial for the interpretation of the neutrino data and for the building of flavour models. 

The purpose of this paper is to study comprehensively the dependence of neutrino 
renormalization effects on all the seesaw parameters, paying special attention to obtaining 
the correct low-energy neutrino measurements. The running of the effective neutrino mass 
matrix below the lightest singlet neutrino mass is generally well under control and large 
renormalization effects can be expected only in the case of degenerate light neutrino masses 
and, in supersymmetric models, for very large values of tan (3 [TT)| ITT] . However, under- 
standing the renormalization effects above and between the heavy neutrino scales fHl E2] is 
much more complicated, since new non-trivial dependences on the heavy neutrino Yukawa 
couplings are introduced. Because the flavour structure of the new contribution to 

1 An up-to-date list of references on this very extensive subject can be found in [T3] . 
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the RGEs can be very different from that due to the effective neutrino mass matrix, large 
effects are possible. Since the couplings Y£> are largely unknown, a typical top-down ap- 
proach taken in previous studies has been to fix the neutrino parameters at Mqut at some 
chosen values, then to run them down to the electroweak scale and demonstrate that, for 
this particular choice, the low-energy neutrino mass matrix is compatible with experimental 
data. 

In this paper we take a bottom-up approach in which we first fix the known low-energy 
neutrino parameters to their measured values, and evaluate renormalization towards higher 
scales consistently in the framework of the minimal supersymmetric seesaw model. In our 
approach, every set of studied neutrino parameters is physical by construction. We pa- 
rameterize the nine high-energy parameters of the seesaw mechanism using the orthogonal 
complex matrix R jlQj, and scan over all the 18 seesaw parameters by generating the un- 
known parameters (including phases) randomly. We run the neutrino parameters up to 
the GUT scale and study the dependence of the renormalization effects on (i) the other 
observable low-energy and (ii) the high-energy parameters. 

We find that significant renormalization effects can occur only when some of the light 
neutrino masses get comparable contributions from two or three heavy neutrinos Nj and/or 
the light neutrino mass scale is at least moderately degenerate. Because the matrix Rij, 
known as the dominance matrix [20] , characterizes the contributions of the heavy neutrino 
Nj to the light neutrino masses z/j, this parametrization turns out to be quite useful for 
the present study. It has been stated in the literature that the solar angle 9 12 usually runs 
more than #13,23- We find that, for light neutrino masses with a strong normal hierarchy, 
exactly the opposite occurs. The quark-lepton complementarity relation [T2*| l2~Tj E21 123] 



turns out to be very stable while, at the same time, the neutrino angles #13 and #23 may 
receive renormalization effects larger than the accuracy of plausible future experimental 
tests. The renormalization of the low-energy oscillation phase 5 is generally enhanced 
compared with the running of mixing angles. Nevertheless, a non-diagonal i?-matrix is 
needed for a large effect also in this case. An interesting feature is the possible crossing of 
light-neutrino mass eigenstates, which is accompanied by discrete changes in the neutrino 
mixing pattern, and is correlated with the i?-matrix parameters. 

Our paper is organized as follows. In Section 2 we present calculational details of our 
study. In Section 3 we present and discuss our results. Finally, some conclusions are drawn 
in Section 4. 

2 Running of Neutrino Parameters in the MSSM 

The superpotential of the minimal supersymmetric standard model (MSSM) with singlet 
(right-handed) heavy neutrinos is given by 



9c + #12 



7T 



(1) 



4 : 




(2) 
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where the Yukawa matrices Y are general complex 3x3 matrices and the 3x3 heavy 
neutrino mass matrix M is symmetric. The Yukawa matrices can be diagonalized by 
bi-unitary transformations Y D = U^YV, where V, U refer to the rotation of the left- and 
right-chiral fields, respectively. In the case of the symmetric matrix M, U = V*. To explain 
the neutrino data naturally, the hierarchy in M should preferably be of the same order as 
the square of the hierarchy in Y v [21]. We therefore assume hierarchical heavy- neutrino 
masses: M 1 < M 2 < M 3 . 

Integrating out all the heavy singlet neutrinos, one gets the usual dimension-5 effective 
operator 

£ = ~kLLH 2 H 2 , (3) 

which after electroweak symmetry breaking gives masses to the light neutrinos: 

m„(/z) = k([i)v 2 sin 2 (3, (4) 

where /i is the renormalization scale, v = 174 GeV and tan/? = v 2 /v\ is the ratio of 
the v.e.v.'s of the corresponding Higgs doublets. Above the heaviest neutrino mass scale, 
> M 3 , the light-neutrino mass matrix reads 

mM = Yj(fi)M~ 1 (fi)Y u (fi)v 2 sin 2 (3. (5) 

Between the heavy-neutrino mass scales, Mi < M 2 < M3, there exist a series of effective 
theories with, in general, n active heavy neutrinos. The tree-level matching conditions 
between these theories at the neutrino thresholds are 

K 9f\M n = Kgf\ Mn + iYy) ng M n (Y V ) ng \ Mn , (6) 

where (n) is the number of heavy neutrinos not integrated out. In general, the light-neutrino 
mass matrix can be written as 



An) (n) m {n) (n) \ o 

Vk + YJM^YA v 2 sin 2 p. (7) 



Since m v and YjY e can be diagonalized with the unitary matrices V u and V e , respectively, 
the mixing matrix observable in the low-energy experiments is 

Vmns = V}V V . (8) 

While m u contains 9 physical parameters, Y v and M together contain 18 parameters. 
The missing 9 parameters crucially affect the physical observables. These include, for ex- 
ample, renormalization-induced lepton-number-violating processes |23 EE EE HI I2Z] and 
electric dipole moments |EJ EH] in the supersymmetric seesaw model, as well as the renor- 
malization of the light-neutrino parameters [T^j. Therefore, to study the dependence of the 
renormalization of eq. ((7J) on Y u , we parametrize Y v with the complex orthogonal matrix 
R JHI- In the basis in which M is diagonal, we write 

(n) (n) , (n) -, 

Y v = (M D )5 R ( m ^)5 Vl (^sin/5)- 1 , (9) 
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where the matrix R is parametrized in terms of three complex angles 9f 2 i 0^ and 



c 12 c 13 6 12 c 13 b 13 

^R pR r?R r?R s^R fiR s^R c?R 



p — I n R Q R C R C R n R n R n R Q R Q R Q R Q R n R I / 1 n\ 

11 — | c 23 6 12 i, 23 i, 13 c 12 c 23 c 12 6 23*13*12 6 23 c 13 I > \ 1K) ) 

„R „R „R Q R „R C .R n R n R C R C R n R n R I 

6 23 6 12 c 23 6 13 c 12 6 23 c 12 L 23*13*12 c 23 c 13 / 

(n) 

where = sin 6^- and = cos 6^. Since Y v and M are renormalized, obviously also R, 
which consists of n rows, runs with energy. The RGEs for Y v and M can be found in [26] . 

(n) 

and the RGEs for R were calculated in [29] . Using these, R has to be evaluated at every 
heavy neutrino threshold when the matching is performed. 

The scale dependence of the effective /combined quantities in (J2J) is characterized by the 
differential equation (SHI UHj 

(n) 

odX j, M>) „(") W , (n) (n) 

levr 2 ^ = (r e tr e + rjnf x + x(yX + rjnf + 

(n) (n) (n) 

(2Tr(Yjy i/ + 3Y^Y U ) _ 6/5^ 2 - 6g 2 2 )X, (11) 

(n) 

where X = k, Yj M~ l Y v . Notice that below the Mi scale Y v = 0. Therefore one expects 
large renormalization effects to occur above the heavy-neutrino thresholds for two reasons. 
First, the Yukawa couplings Y v can be large. Secondly, the flavour structure of YJY V can 
be very different from the flavour structure of YJY e and k. Both those effects can be traced 
back to the values of R via eq. (jHJ). Approximate analytical solutions to eq. (fTTjl have been 
derived in ^3], which allow one to understand the generic behaviour of the renormalization 
effects. However, due to enhancement /suppression factors and possible cancellations, the 
exact numerical solutions may differ considerably from those estimates. 



3 Results for Normally- Ordered Light Neutrinos 

In this Section we present the results of our study for the case of normally-ordered light- 
neutrino masses, using the following strategy. We start at Mz, where we fix the measured 
light-neutrino parameters as Am 2 2 = 8. x 10~ 5 eV 2 , Am^ = 2.2 x lCT 3 eV 2 , tan 2 12 = 0.41, 
sin 6 2 s = 0.7 and sin #13 = 0.05. We then generate randomly the lightest neutrino mass, the 
heavy neutrino masses, all the low-energy phases and the initial values for the parameters 
in the R matrix. We run all the relevant quantities up to Mqvt using the 1-loop RGEs 
for the minimal supersymmetric seesaw model |2SHSH|- At every heavy- neutrino threshold 

(n) 

we perform the tree-level matching according to eq. (jBJ). To calculate the values of Y v we 
use the renormalized values of the R and M. At Mqut we calculate the renormalized 
light-neutrino parameters. We always keep the ordering of the light neutrino masses fixed, 
nil < m 2 < m 3 for the normal and m 3 < mi < m 2 for the inverse ordering. Because of 
that the physical range for 612 extends up to tt/2. 

In order to accentuate the renormalization effects due to the low-energy neutrino pa- 
rameters and the parameter matrix R, we do not consider degenerate light neutrinos and 
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we assume an upper limit m\ < 0.1 eV on the lightest neutrino mass. Although the present 
most stringent limit on the overall light-neutrino mass scale scale coming from astrophysics 
and cosmology jSTj is a factor of 2 to 3 weaker, such precision can easily be achieved in the 
future cosmological experiments [Zj. We also minimize the renormalization effects of large 
tan/3 studied in by working with the relatively small value of tan/3 = 5. Instead, we 
study how the large values of Y v affect the renormalization effects. 

3.1 Renormalization of the Mixing Angles 

We start by studying the running of the light-neutrino mixing angles. In Fig. ^ we plot the 
neutrino mixing angles % at Mqut as functions of the lightest neutrino mass m,\{Mz) for 
R(Mz) = 1 (left panel) and for randomly generated R (right panel). In all the figures the 
neutrino mass parameters are presented in units of eV. The mixing angles #12, #13 and 9 2 3 are 
represented by green (light), blue (dark) and red (medium) dots, respectively. For R = 1 
the mixing angles practically do not run: only 8 12 may change a little for light-neutrino 
masses close to m\ = 0.1 eV. This is because, for moderate degeneracy, the renormalization 
of Q\2 is enhanced by a factor m\ / Am\ 2 . This effect would have been larger for larger values 
of tan (3 and mi. 

Turning to the results for the randomly-generated values of R, a certain pattern of 
renormalization effects emerges. If mi > ^ Am\ 2) very large changes in the mixing angles 
may occur. Although # 12 tends to change most due to the above-mentioned enhancement 
factor, also #13 and 823 may gain almost any value. We see that the examples of extreme 
running considered in the literature are due to having at least a moderately degenerate 
mass spectrum. 

An interesting set of points in Fig. ^ are those gathered around 612 ~ 60° in the region 
mi > a/ Am Those correspond to the level crossing of two light-neutrino mass eigenvalues 
mi and UI2 due to renormalization. Because by definition mi < m 2 , this causes a discrete 
jump in the value of #12. In the standard parametrization of VyiNSi and with small #13, this 
implies sin#i 2 — > cos#' 12 and, consequently, 9[ 2 = 90° — ^12- As seen i n Fig. EJ this effect is 
smeared by strong running of 9i2 and also #13. 

In contrast to the previous discussion, if the mass spectrum is strongly hierarchical: 
mi < \J Amf 2 , the solar angle is much more stable than the mixing angles #13 and #23- The 
latter may vary through a range of almost 10°, which is more than the expected precision 
of future experiments. Moreover, the widths of the #13 and #23 bands in Fig. ^ do not de- 
pend on the initial values of the angles Thus, discrimination between different flavour 
models may be possible in principle in the future, if one takes into account renormalization 
effects. We also note that, for hierarchical light-neutrino masses, the quark- lepton comple- 
mentarity relation 9c + #12 = would be maintained with high accuracy at every scale, 
independently of the unobservable neutrino parameters. 

We now study the origins of the effects due to R. In the left panels of Figs. 121 E]and|I] 
we plot the distributions of the neutrino mixing angles 9ij(McuT) as functions of mi(Mz) 
for randomly generated complex parameters 9^, 9f z and 9^, respectively. In each figure 
the other parameters in R are set to zero. The same neutrino mixing angles are plotted in 
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Figure 1: Neutrino mixing angles at Mqut cls functions of the lightest neutrino mass at 
mi(Mz) for R(Mz) = 1 (left panel) and randomly generated R (right panel). 



the right panels of Figs. 121 El HI as functions of the absolute values of the corresponding R 
matrix parameters |#^|, |#{|| and |#H|, respectively. 

We see in Fig. [2] that non-zero values of 9± 2 affect mostly the renormalization of #12. 
A large running effect requires also that the overall light neutrino mass scale be high. On 
the other hand, non-zero values of 8± 3 affect mostly the running of #13 and 823, as seen in 
Fig. El Again, a relatively high overall light-neutrino mass scale is required for a significant 
effect. We also learn from Fig.EJthat the level crossing of light mass eigenvalues is induced 
by non-zero 9± 3 , which strongly affects the running of mi. The parameter 9 23 affects only 
the running of #13 and #23 . Fig. 0] shows an interesting feature: in this case the running of 
813 and 623 does not depend on mi, and significant renormalization effects can be obtained 
also for very hierarchical light neutrinos. 

Comparison of the left and right panels in Figs. 121 El and 0] reveals how the renormaliza- 
tion effects depend on the magnitude of the particular parameter 9^. Interestingly, in all 
the cases the dominant running occurs in the region \6^\ ~ 0(1). We recall that R can be 
interpreted to be a dominance matrix [20J, i.e., it shows which heavy neutrino contribution 
dominates in the mass of a particular light neutrino. Therefore our results imply that, in 
order to have significant renormalization effects, at least two heavy neutrinos must con- 
tribute to one particular light neutrino mass in approximately equal amounts. If the light 
neutrino masses are dominated by one heavy neutrino each, no large running is possible 
unless the light neutrinos are degenerate in mass and tan (3 is large. 



3.2 Renormalization of Masses 

The observed hierarchy in the light-neutrino masses, a/ Amf 2 / Am^ = 0.18, is milder than 
expected in many flavour models. In GUTs with the simplest scalar sector, the neutrino 
Yukawa couplings are equal to the up-quark Yukawa couplings at M GUT . Contrary to that, 
phenomenology at the low scale seems to indicate that the neutrino hierarchy is more similar 
to the less hierarchical down-sector Yukawa couplings, rather than those in the up sector. 
However, at the moment the lightest neutrino mass and the hierarchy in the heavy-singlet 
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Figure 2: Neutrino mixing angles at Mqut as functions of the lightest neutrino mass 
mi(Mz) for random 6f 2 fl e ft panel), and as functions of \9^\ (right panel). The rest 
of the parameters in RiMz) vanish. 
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Figure 3: Neutrino mixing angles at Mqut as functions of the lightest neutrino mass 
mi(Mz) for random 9± 3 (left panel), and as functions of \9± 3 \ (right panel). The rest 
of the parameters in RiMz) vanish. 
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Figure 4: Neutrino mixing angles at Mqut as functions of the lightest neutrino mass 
mi(M z ) for random 9 23 (left panel), and as functions of \9^\ (right panel). The rest 
of the parameters in R(Mz) vanish. 
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Figure 5: The distribution of a/ Am^ versus 
and for randomly generated R (right panel). 



\P±m\ 2 at Mqut for R{M Z ) = 1 (left panel) 



sector are unknown, introducing large uncertainties into such considerations. Therefore it 
is interesting to study also how the masses of the light neutrinos evolve with energy. 

In Fig. we plot the distributions of a/ Amf 2 and a/ Am^ at M GUT for R(M Z ) = 1 
(left panel) and for randomly generated R (right panel). The red dot denotes the starting 
point at Mz from which value all the other points are generated. The hierarchy at the 
GUT scale tends to be larger than at low energies, although the opposite is also possible for 
a few points. While for R = 1 both mass differences tend to increase, for random R they 
may also decrease. The abundant points with smaller values of a/ Amf 2 in the right-hand 
plot correspond to non-zero values of 9^ 3 . This parameter affects the Yukawa couplings of 
first generation in such a way that the (12) mass difference may run considerably. 

3.3 Renormalization of CP Observables 

Our approach in this study is to fix the known neutrino parameters and to vary the un- 
known ones randomly. At the moment, the only CP- violating observable in the neutrino 
sector what we have information about is the baryon asymmetry of the Universe, assuming 
that the cosmological baryon asymmetry is generated via leptogenesis [321 • In this case, it 
is possible to constrain one combination of the 6 CP-violating phases in the neutrino sec- 
tor. However, because one can vary the remaining 5 combinations (and also the unknown 
CP-conserving neutrino parameters), one cannot make any firm prediction for the neutrino 
parameters [2Z1- As the predictions for other renormalization- induced CP- violating observ- 
ables such as the electric dipole moments of the charged leptons are orders of magnitude 
smaller than the present experimental bound ^JEHj, no firm constraints come from this 
sector either. Although Ref. |EJ argues that some systematic renormalization effects in 
leptogenesis are possible due to the running of the effective light neutrino mass matrix, 
those are already taken into account in the systematic study of 

At the moment, the most realistic possibility seems to be that of measuring the MNS 
phase 5 in future oscillation experiments. The value of this phase, however, is presently 
unknown. In the left plot of Fig. we present the values of 5{Mqut) against the initial 
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Figure 6: The MNS phase at high energy S(Mgut) a function of its low-energy value 
5(Mz) for R(Mz) = 1 (left panel), and 5{Mgut) — 8(Mz) as a function of mi for randomly 
generated R (right panel). 



values of the phase for R(Mz) = 1. As seen in the figure there is practically no running of 6 
in this case. The situation changes considerably for randomly generated i?-matrices, as seen 
in the right panel of Fig. |H1 where we plot the change of the phase, 5(Mgut) — S(Mz), as a 
function of the lightest neutrino mass m\. The running of the MNS phase can be numerically 
significant and, apart from high values of m\, almost independent of the lightest neutrino 
mass. This behaviour resembles the running of 613^3 in Fig. ^ and can be traced back to 
the l/#i3 enhancement of the running of 5 ^3] 2 ■ 



4 Renormalization Effects for an Inverted Hierarchy 
of Light Neutrino Masses 

We have repeated the earlier analyses also for an inverted hierarchy of light neutrino masses. 
Because the solar-neutrino mass scale is now higher than the atmospheric one, the Yukawa 
couplings of the first generations are generally larger. Therefore, all the effects related 
to the renormalization of the solar parameters are generally enhanced. This can be seen 
in Fig. [7J where we plot the neutrino mixing angles at Mqut as functions of the lightest 
neutrino mass rri3(Mz) for R(Mz) = 1 and for randomly-generated R. If R(Mz) = 1, the 
mixing practically does not run even in the inverted-hierarchy case. The exception is in 
the high- mass region rri3(Mz) ~ 0.1 eV, when the mass eigenvalues cross and the moderate 
degeneracy causes discrete jumps in the mixing angles. As expected, for the random choice 
of R, the angle (9 12 runs very strongly. The level-crossing stripe around 812 ~ 60° exists for 
all values of 7713, while the angles #13 and 623 run only for large values of 7713. 

2 The points around ±2ir in Fig. actually correspond to small running modulo 2tt. 
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Figure 7: Neutrino mixing angles at Mqut as functions of the lightest neutrino mass 
m^{Mz) for R(Mz) = 1 (left panel) and for randomly generated R (right panel). 



5 Discussion and Conclusions 

We have studied how the RGE running of neutrino mixing depends on the unknown seesaw 
parameters. We have taken a bottom-up approach in which we fix the known low-energy 
neutrino parameters to their measured values. Parametrizing the nine high-energy param- 
eters of the seesaw mechanism via the dominance matrix R, we have scanned over all the 
18 seesaw parameters by generating the unknown parameters randomly. The fact that the 
matrix Rij measures the heavy neutrino Nj contribution to the light neutrino z/j means 
that this parametrization is particularly valuable for this study. We have compared the 
results for random R elements with the simple case R = 1. 

We have found that significant running effects can occur only when some of the light- 
neutrino masses have comparable contributions from more than one heavy neutrino iVj, and 
the light-neutrino mass scale is at least moderately degenerate. For a normal hierarchy of 
neutrino masses, the complementarity relation (1) between neutrino and quark mixing 
angles evolves very little between the GUT scale and the electroweak scale. However, 
the other oscillation angles #13^3 run rather more than 6*12 and also beyond the expected 
measurement errors. In certain cases, we observe level crossing in the light-neutrino mass 
eigenstates that is accompanied by jumps in the oscillation angles. The running of the 
CP- violating oscillation phase S is strong for random R but insignificant for R = 1. 

The analysis presented here complements the top-down approach often adopted else- 
where. It reveals some of the pitfalls in inferring properties of the neutrino mass matrix 
generated at the GUT scale from current low-energy measurements alone, in the absence 
of supplementary theoretical or phenomenological input. We hope that these results may 
serve as useful aids in the attempt to understand the neutrino mass matrix, which has 
already revealed several surprises. The results presented here demonstrate that our low- 
energy measurements are far from telling us the whole story. 
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